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i. 


ABSTRACT 


The  effects  of  topographic  features  on  Rayleigh  wave  propagation  and 
scattering  are  investigated  in  the  laboratory  using  three-dimensional  ultrasonic 
models.  Starting  from  simple  steps,  different  topographic  features  are  modeled. 
The  effects  of  these  features  on  Rayleigh  wave  transmission  and  scattering  are 
examined  as  a  function  of  wavelength  and  as  a  function  of  angle  of  incidence. 

In  general,  backscattered  or  reflected  Rayleigh  waves  are  small  compared  to 
transmitted  waves.  A  significant  fraction  of  the  Rayleigh  wave  energy  is 
scattered  into  body  waves.  To  complement  the  laboratory  studies,  two- 
dimensional  finite  difference  calculations  are  carried  out.  These  show  the 
details  of  the  scattering  of  incident  Rayleigh  waves  at  a  step  discontinuity  into 
P,  S,  and  Rayleigh  waves.  More  energy  scatters  into  S  waves  than  into  P  waves. 

Using  three-dimensional  laboratory  models  of  steps,  the  dependence  of 
wave  scattering  on  angle  of  incidence  is  studied.  Transmission  and  reflection 
coefficients  (transmitted  or  reflected  energy/incident  energy)  computed  from 
spectral  ratios  vary  strongly  with  incidence  angle.  At  wavelengths  equal  to 
twice  the  step  height,  the  fraction  of  incident  energy  scattered  into  body  waves 
ranges  from  more  than  90%  at  normal  incidence  to  about  zero  at  near-grazing 
incidence.  At  each  angle,  transmission  coefficients  vary  strongly  with 
frequency.  Because  of  frequency-dependent  phase  shifts,  the  transmitted  and 
reflected  waves  are  distorted. 

The  effect  of  the  steps  on  the  propagation  of  Rayleigh  waves  is 
demonstrated  by  convolving  synthetic  dispersed  wave  trains  with  the  impulse 
response  of  the  scale  models.  The  ocean-continent  margin  of  the  western 
United  States  is  modeled  as  a  60°  ramp  scaled  to  60  km  height.  The  Tibetan 
Plateau  is  modeled  as  a  broad  mesa  scaled  to  40  km  height.  In  both  models  the 
azimuthal  dependence  of  transmitted  Rayleigh  waves  is  similar  to  that  observed 
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1  INTRODUCTION 


Lateral  crustal  heterogeneities  have  an  important  influence  on  the 
propagation  of  surface  waves.  In  Figure  1,  long-period  surface  waves  recorded 
from  an  isotropic  source  (mb  =  7  explosion)  in  Eastern  Kazakhstan  are  shown 
at  Eurasian  stations.  Surface  waves  crossing  Tibet  and  the  Himalayas  are 
severely  attenuated.  These  observations  are  similar  to  those  of  Bird  and  Toksoz 
(1977).  Similarly,  Rayleigh  waves  recorded  from  a  nuclear  explosion  at  the 
Nevada  Test  Site  are  strongly  attenuated  as  they  cross  the  Pacific  margin  as 
shown  in  Figure  2.  In  recent  years  studies  have  concentrated  on  regional 
variations  of  the  lithosphere  and  mantle  structures.  The  effects  of  such 
heterogeneities  on  surface  waves  have  been  considered  primarily  in  terms  of 
deviation  of  propagation  paths  from  the  "great  circle"  path  and  multipathing 
effects  (Yomogida  and  Aki,  19B5).  Yet  as  Figures  1  and  2  show,  major 
topographic  features  and  sharp  structural  contrasts  at  ocean-continent 
boundaries,  active  subduction  zones  and  continental  collision  belts  could  have 
significant  effects  on  surface  wave  amplitude  attenuation  due  to  scattering. 

Higher  modes  of  surface  waves  are  affected  by  propagation  across  lateral 
heterogeneities  as  well.  Rusaikin  et  al.  (1977)  showed  that  lg  was  absent  from 
paths  crossing  Tibet,  and  proposed  three  possible  explanations  for  this 
absence:  (l)  abrupt  changes  in  crustal  structure  at  the  plateau  margins,  (2) 
unusual  velocity  structure  beneath  the  plateau,  and  (3)  low-#  regions  within 
the  plateau.  Chinn  et  al.  (1980)  studied  the  propagation  of  Lg  in  the  Andes,  and 
suggested  that  the  primary  factor  in  the  efficient  propagation  of  lg  is  the 
orientation  of  the  path  relative  to  the  structural  trend  of  the  mountains: 
efficient  lg  propagation  parallel  to  the  structural  trend,  inefficient  propagation 
perpendicular  to  that  trend.  Baumgardt  (1985)  found  a  similar  behavior  for 
propagation  across  the  Urals,  The  purpose  of  our  study  is  to  understand  the 
surface  wave  scattering  across  topographical  and  structural  heterogeneities 


using  ultrasonic  modeling.  In  ultrasonic  laboratory  studies  we  observe  exactly 
the  same  trend  of  waves  attenuating  strongly  across  structural  heterogeneities. 
Figure  3  shows  laboratory  results  of  Rayleigh  wave  propagation  in  a  model  of 
Nevada  Basin  and  Range  structure.  Note  that  from  a  point  source  at  S,  the 
Rayleigh  wave  propagating  parallel  to  the  structural  trend  (R \)  is  much  larger 
than  those  that  go  across  the  trend  ( R%  and  /?g).  In  the  direction  parallel  to  the 
trend,  the  wave  shows  relatively  little  scattering,  while  normal  to  the  trend  wave 
shapes  have  been  altered  and  attenuated  because  of  strong  scattering.  These 
effects  are  investigated  in  detail  in  the  next  sections  using  both  two-  and 
three-dimensional  models. 
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FIGURE  CAPTIONS 


Figure  1.  Long  period  vertical  component  surface  waves  recorded  at  WWSSN 
stations  (triangles)  from  Mb  =  7  explosion  of  11/30/69  (Eastern  Kazakhstan). 

Figure  2.  Rayleigh  waves  crossing  the  continent-ocean  boundary;  seismograms 
from  the  nuclear  explosion  Boxcar,  recorded  at  Berkeley  (BRK)  and  Ocean 
Bottom  Seismometer  (OBS)  in  the  Pacific.  At  top  locations  of  the  Nevada  Test 
Site,  BRK  and  OBS  stations  and  the  path  are  shown.  Traces  (a)  and  (b)  are  raw 
(unnormalized)  seismograms  at  BRK  and  OBS,  respectively.  The  two  traces  at 
(c)  are  filtered  and  gain-normalized.  The  solid  trace  is  BRK  and  the  dashed 
trace  OBS.  Note  that  amplitude  decreases  by  about  a  factor  of  3  when  the  wave 
crosses  the  continent-ocean  margin  (after  McGarr,  1969). 

Figure  3.  Data  taken  with  valleys  filled  with  two  parallel  layers.  Effects  of 
topography  and  layered  valley  on  Rayleigh  wave  scattering  and  dispersion. 
Model  of  Dry  Lake  Valley,  Nevada  (scale  1:250,000). 
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H.  TWO-DIMENSIONAL  ULTRASONIC  EXPERIMENTS  AND  FINITE 
DIFFERENCE  CALCULATIONS 

11.1.  LABORATORY  STUDIES  OF  SURFACE  SCATTERING 

Two  and  three-dimensional  ultrasonic  laboratory  experiments  have  been 
carried  out  on  the  scattering  of  Rayleigh  waves  from  surface  topographic  and 
structural  heterogeneities.  Two  dimensional  experiments  were  done  using 
nickel  and  aluminum  plates  with  topographic  discontinuities  modeled  as  steps, 
ramps,  rectangular  and  triangular  shaped  "mountains”  (Nathman,  1980;  Toksoz, 
1983).  These  experiments  extended  the  results  of  previous  studies  (de 
Bremaecker,  1958;  Kato  and  Takagi,  1956;  Knopoff  and  Gangi,  1960;  Lewis  and 
Dally,  1970;  Martel  et  at.,  1977;  Pilant  et  a l.,  1964).  Figures  4,  5  and  6  show 
examples  of  incident,  reflected  (back-scattered)  and  transmitted  (forward- 
scattered)  Rayleigh  waves  across  up  and  downsteps  and  "mountains."  Note  that 
transmitted  and  reflected  waves  account  for  a  small  fraction  of  incident 
Rayleigh  wave  energy.  A  large  fraction  of  Rayleigh  wave  energy  is  scattered  into 
body  waves.  Figures  7a,b  show  the  nature  of  body  waves  generated  due  to 
Rayleigh  wave  scattering  at  a  step.  Note  that  more  energy  is  scattered  into  S 
waves  and  both  P  and  S  scattering  patterns  are  complex. 

II.2  FINITE  DIFFERENCE  MODELING  OF  SURFACE  WAVE  SCATTERING 

To  investigate  details  of  the  scattering  process,  two-dimensional  finite 
difference  calculations  are  carried  out  for  Rayleigh  waves  incident  on 
topographic  steps  (Toksoz,  1983;  Fuyuki  and  Nakano,  1984;  Stephen  et  a l,, 

1985).  In  a  perfectly  elastic,  two-dimensional,  homogeneous,  isotropic  medium, 
the  wave  propagation  for  displacement  u  can  be  written 


U  ti  -A ,  u  «  +AZ  u  „  +AS  u  „ 


(1 


where  the  displacement  vector  u  equals 

u  =  (  u(x,z,t)  ,  w(x,z,t )  ) 

and 
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where  a  and  /?  are  P  and  S  wave  velocities.  The  boundary  conditions  for  the 
free  surface  are  the  vanishing  normal  and  shear  stresses.  These  conditions  can 
be  represented  in  the  case  of  an  horizontal  free  surface  by 


u*+B1ux=  0  (2 

and  in  the  case  of  a  vertical  free  surface  by 


where 
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The  finite  difference  method  is  used  to  solve  equation  (1)  subject  to  the 
boundary  conditions  (2)  and  (3)  and  to  the  initial  conditions  which  will  be 
described  later.  The  procedure  followed  is  similar  to  that  used  by  Fuyuki  and 


Matsumoto  (1980),  and  Munasinghe  and  Farnell  (1973). 


For  the  case  of  propagation  in  an  infinite  homogeneous  media,  the  explicit 
finite  difference  formulation  is  stable  only  if: 


Lt<. 


min(Az.Az) 

vTot2+02) 


(4) 


when  Az  and  Az  are  the  grid  dimensions  in  the  horizontal  and  vertical 
directions  and  A t  is  the  timestep.  The  principal  cause  of  inaccuracy  in  finite 
difference  calculations  in  slowly  varying  media  is  grid  dispersion.  If  the  grid 
increments  (Az,  Az)  are  too  large,  low  frequencies  will  travel  faster  across  the 
grid  than  high  frequencies,  causing  apparent  dispersion.  This  result  is  generally 
true  for  compressional  waves.  For  shear  waves  the  dispersion  relation  is  more 
complex.  For  some  values  of  Poisson's  ratio  and  propagation  direction  high 
frequencies  may  travel  faster  than  low  frequencies  (Bamberger  el  al.,  1980; 
Trefethen,  1982;  Sato  and  Jshihara,  1983).  Care  was  taken  to  assure  that  both 
of  the  potential  problems  were  correctly  handled. 


In  our  calculations  we  generate  the  initial  Rayleigh  wave  with  its 
eigenvelocity  Vr.  The  incident  waveform  is  chosen  so  that  the  vertical 
displacement  to  on  a  free  surface  has  the  shape  of  a  Ricker  wavelet  (Ricker, 
1953).  The  waveform  in  the  wavenumber  spectrum  is 


S(k  )=(fc  /  ko)2exp[  1  —(k  /  k  0)2]  (5) 

where  denotes  the  wavenumber  corresponding  to  the  maximum  amplitude, 

and  the  dominant  wavelength  Xo  is  defined  by  Xo=2  T~-  This  waveform  has  been 

k0 

used  before  by  Munasinghe  and  Farnell  (1973),  and  Fuyuki  and  Matsumoto 
(1980).  To  eliminate  the  edge  effects,  absorbing  boundaries  were  placed  at  the 
sides  and  the  bottom  of  the  model  following  the  formulation  of  Clayton  and 
Engquist  (1977)  and  Reynolds  (1978). 
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In  order  to  simulate  laboratory  models,  the  finite  difference  calculations 
are  carried  out  using  the  model  parameters  with  P  wave  velocity  6  km/sec,  S 
wave  velocity  3.46  km/sec  and  Rayleigh  wave  velocity  3.18  km/sec.  This  case 
corresponds  to  a  Poisson  ratio  of  0.25.  The  center  frequency  of  the  wavelet  is 
0.318  MHz  and  the  corresponding  wavelength  is  10  mm.  In  numerical  modeling 
we  use  25  grid  points  per  wavelength,  which  gives  Ax=Az  =0.4  mm  The  At  value 
is  0.047  (j,  sec. 

The  ultrasonic  models  have  corners  at  90°  and  270°  angles.  A  partially  off- 
centered  finite  difference  scheme  has  been  employed  by  other  authors 
(Alterman  and  Rotenberg,  1969;  Munasinghe  and  Farnell,  1973;  Fuyuki  and 
Matsumoto,  1980)  to  deal  with  such  corners.  The  first  two  methods  give  an 
accuracy  of  the  order  of  Ax,  and  the  Fuyuki  and  Matsumoto  method  an 
accuracy  on  the  order  of  (Ax)2.  Fuyuki  and  Matsumoto  (1980)  and  Fuyuki  and 
Nakano  (1984)  show  that  the  relative  error  in  finite  difference  calculations  is  ± 

1  per  cent  for  the  horizontal  and  vertical  displacements  located  close  to  the 
free  surface.  Below  one  wavelength  the  relative  errors  increase  to  about  ±4  per 
cent.  Since  we  follow  the  same  approach  as  Fuyuki  and  Matsumoto,  the 
accuracy  of  our  calculations  should  be  similar  to  the  above  values. 

In  order  to  separate  the  scattered  P  and  S  waves,  the  divergence  (V.u)  and 
curl  (V>§i)  of  the  displacement  field  is  calculated  at  each  point  at  designated 
times.  Figures  8  and  9  show  snapshots  of  the  motion  at  different  times  for 
Rayleigh  waves  incident  on  a  downstep  and  upstep  surface  discontinuity.  Note 
that  the  dilatation  component  of  the  scattered  field  is  much  smaller  than  the 
rotational.  The  complex  patterns  of  scattering  of  Rayleigh  to  P  and  S,  and  in 
turn  the  scattering  of  body  waves  are  clearly  shown  in  these  color  figures. 

These  results  show  in  greater  detail  what  was  observed  in  the  laboratory 
experiments  (Figures  4,  5,  6,  7a,  and  7b). 
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A  case  of  special  interest  is  the  Rayleigh  wave  scattering  from  a 
rectangular  surface  feature  (i.e.,  a  mesa).  The  experimental  results  showed 
that  two  Rayleigh  waves  are  generated  by  the  mesa:  the  first  by  the  Rayleigh- 
shear-Rayleigh  scattering,  the  second  by  Rayleigh  wave  transmission  following 
the  sfrface  contour  (Figure  6,  especially  second  trace  from  bottom). 

Finite  difference  calculations  for  this  case  are  shown  in  Figures  10  and  11. 
Figure  10  shows  the  particle  displacement  vectors.  Frames  c,  d,  e,  and  f  show 
how  Rayleigh-shear-Rayleigh  scattering  takes  place.  The  complexity  of 
displacements  due  to  interference  in  the  mesa  is  clearly  illustrated.  In  Figure 
11.  the  curl  of  displacement  (shear  strain)  on  the  right  hand  side  demonstrates 
a  case  of  Rayleigh-shear-Rayleigh  scattering.  Rayleigh  waves  "transmitted"  by 
this  mechanism  are  larger  than  Rayleigh  waves  that  follow  the  surface. 


In  summary,  these  two-dimensional  results  demonstrate  that  high 
frequency  surface  wave  propagation  across  topographic  irregularities  can  be 
complex.  Multiple  scattering  between  Rayleigh  and  body  waves  can  take  many 
forms  based  on  topography  and  wavelength.  In  the  next  section  we  show  the 
results  from  simple  three-dimensional  models. 
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FIGURE  CAPTIONS 


Figure  4.  Seismograms  of  reflected  and  transmitted  two-dimensional  Rayleigh 
wave  pulses  incident  on  an  upstep. 

Figure  5.  Seismograms  of  reflect  "d  and  transmitted  two-dimensional  Rayleigh 
waves  incident  on  a  downstep. 

Figure  6.  Seismograms  of  reflected  and  transmitted  pulses  for  rectangular 
"mountains"  with  different  geometries. 

Figure  7a.  Receiver  positions  used  to  measure  azimuthal  scattering  of  step 
(two-dimensional). 

Figure  7b.  Seismograms  taken  at  receiver  positions  as  shown  on  Figure  7a. 

Figure  8a,b.c.  Time  sequence  §£  divergence  (left  side)  and  curl  (right  side)  of 
displacement  for  Rayleigh  wave  incident  on  a  downstep. 

Figure  9.  As  for  Figure  8.  but  for  upstep. 

Figure  lOa-i.  Time  sequence  of  displacement  vectors  depicting  the  scattering 
of  a  surface  wave  for  the  rectangular  "mountain"  case.  See  text  for  discussion 

Figure  lla.b.  As  for  Figure  8,  but  for  a  rectangular  "mountain". 
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ID.  THREE  DIMENSIONAL  EXPERIMENTS 


III.1  INTRODUCTION 

Rayleigh  wave  propagation  in  strongly- varying  structures  is  a  complex 
process  well  suited  to  analysis  by  ultrasonic  modeling  methods.  Until  recently, 
ultrasonic  modeling  of  surface  wave  transmission  and  scattering  has  been 
limited  to  two  dimensions.  In  this  chapter  we  report  on  a  series  of  ultrasonic 
modeling  experiments  investigating  the  propagation  and  scattering  of  Rayleigh 
waves  by  simple  three-dimensional  topographic  features. 

Rayleigh  waves  propagating  in  two  dimensions  through  simple  steps  or 
around  corners  and  in  wedges  have  been  studied  experimentally  by  a  number  of 
investigators  (de  Bremaecker,  1958;  Kato  and  Takagi,  1956;  Lewis  and  Dally, 

1970;  Knopoff  and  Gangi,  1960;  Pilam  et  al.f  1964;  Martel  et  at.,  1977;  and 
Nathman,  1980).  Detailed  'snapshots'  of  the  complete  stress  field  of  Rayleigh 
waves  scattered  in  wedges  have  been  obtained  using  photoelastic  methods 
(Lewis  and  Dally,  1970).  Chapter  II  of  this  report  presents  our  results. 

Development  of  theoretical  and  numerical  methods  of  synthesizing 
Rayleigh  waves  scattered  from  simple  structures  has  typically  lagged  behind 
experimental  results.  Rayleigh  wave  transmission  and  reflection  coefficients 
across  vertical  boundaries  in  two  dimensions  have  been  calculated  using 
approximate  variational  methods  (e.g.  McGarr  and  Alsop,  1967).  Recently, 
synthesis  of  complete  seismograms  of  Rayleigh  waves  reflected  from,  and 
transmitted  through,  a  vertical  step  has  been  performed  using  the  finite- 
difference  method  (Toksoz  1983;  Fuyuki  and  Nakano,  1984;  Chapter  II).  In  three 
dimensions,  numerical  computations  are  limited.  Approximate  methods  to 
estimate  transmission  and  reflection  coefficients  of  Rayleigh  waves  incident  at 
steep  angles  on  a  vertical  boundary  have,  however,  been  successfully  applied 
(Malichewsky,  1976;  Chen  and  Alsop,  1979). 


Scattering  of  Rayleigh  waves  in  three  dimensions  can  be  investigated  in  the 
laboratory  using  ultrasonic  methods.  Rayleigh  wave  propagation  in  three- 
dimensions  has  already  been  studied  using  realistic  three-dimensional  scale 
models  of  surface  topography  (e.g.  Toksoz,  1983).  In  such  models,  however, 
even  a  simple  input  pulse  is  severely  distorted  by  the  complex  structure. 

Subtle  effects  of  topography  on  Rayleigh  wave  propagation  and  scattering  are 
thus  easily  obscured  by  the  complexity  of  the  model.  An  important  step  toward 
an  understanding  of  surface  wave  scattering  in  three-dimensions  is  the 
investigation  of  propagation  in  simple  models  of  canonical  form.  In  this 
chapter,  we  investigate  Rayleigh  wave  scattering  from  step  and  ramp-type 
topographic  discontinuities  and  study  their  implications  for  surface  wave 
propagation  across  major  features  such  as  ocean-continent  boundaries  and  the 
Tibetan  Plateau. 

III. 2  EXPERIMENTAL  METHOD 

The  model  used  in  this  study  consists  of  a  solid  block  (200x200x100  mm^) 
of  aluminum  Steps  of  various  height  and  geometry  are  milled  into  one  of  the 
large  faces  (Figure  12). 

The  ultrasonic  source  and  receiver  each  consist  of  a  2.25  MHz  P  wave 
transducer  coupled  to  a  lucite  wedge.  P  waves  generated  by  the  transducer 
propagate  through  the  wedge  and  excite  Rayleigh  waves  in  the  aluminum  with 
high  efficiency.  Similarly,  the  transducer-wedge  combination  is  an  efficient 
Rayleigh  wave  sensor.  The  radiation  pattern  of  the  transducer-wedge 
combination  is  highly  directional;  at  an  angle  of  10°  from  the  axis  of  the  central 
lobe,  peak  energy  is  reduced  to  less  than  25%  of  the  axial  peak  energy.  The 
combination  of  these  features  makes  the  transducer-wedge  combination 
especially  useful  when  working  with  a  low-loss  medium  such  as  aluminum; 
attenuation  can  be  neglected  in  the  analysis  of  Rayleigh  wave  propagation,  and 


the  observed  Rayleigh  waves  are  uncontaminated  by  body  waves  multiply 
reflected  at  the  edges  of  the  model. 
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Peak  energy  of  Rayleigh  waves  traveling  on  the  aluminum  block  ( Vp  =  6.4 
km/sec;  ys  =  3.2  km/sec;  cr  =  3.0  km/sec)  lies  in  the  range  0.2— 2.0  MHz 
(A=  1.5—15  mm).  In  this  frequency  band  the  spectrum  is  repeatably  obtained 
from  a  given  model  and  source-receiver  configuration.  We  therefore  take  this 
as  the  usable  bandwidth  of  the  modeling  experiments.  Because  the  model  is 
essentially  a  half-space,  only  fundamental-mode  Rayleigh  waves  are  observed. 
Since  a  half-space  is  non-dispersive  to  Rayleigh  waves,  the  half-space  Rayleigh 


wave  represents  the  impulse  response  of  the  transducers  and  recording  system. 

When  a  Rayleigh  wave  encounters  a  step  some  of  the  incident  energy  is 
transmitted  through  the  step,  some  is  reflected  from  the  step,  and  some  is 
converted  to  body  waves  (de  Bremaecker,  1958;  Knopoff  and  Gangi,  1960;  Martel 
et  a Z  ,  1977;  Toksoz,  1983).  In  general  the  efficiency  of  transmission  through  or 
reflection  from  a  step  depends  on  the  wavelength  relative  to  the  step  height 
(e.g.  McGarr  and  Alsop,  1967).  In  order  to  facilitate  comparison  of  models  with 
different  step  heights,  therefore,  wavelengths  are  normalized  to  the  step  height. 
The  "normalized  frequency",  /  ,  is  obtained  from  the  scaling  relation 


?  -  -  T  (6) 

where  f  is  the  frequency,  c  is  the  phase  velocity,  A  is  the  wavelength,  and  h  is 
the  step  height.  For  a  given  step  height,  the  range  of  usable  normalized 
frequencies  is  limited  by  the  usable  bandwidth  of  the  modeling  experiments. 

This  is  an  important  consideration  when  the  model  seismograms  are  scaled  to 
correspond  to  real  Earth  seismograms.  Transmission  and  reflection  coefficients 
calculated  from  Rayleigh  waves  traveling  across  steps  of  various  heights 
indicate  that  the  model  scaling  relation  (Equation  6)  is  generally  appropriate 
over  a  wide  range  of  frequencies  and  step  heights. 


We  define  the  the  Rayleigh  wave  energy  transmission  coefficient,  Ef,  as 


M /)  =  7777  .  (7a) 

AoV  ) 

where  Ar(f)  is  the  Fourier  amplitude  spectrum  of  the  transmitted  Rayleigh 
wave,  and  Aq(/)  is  the  amplitude  spectrum  of  the  Rayleigh  wave  propagating 
across  the  aluminum  half-space.  The  reflection  coefficient,  Er,  is  defined  in  an 
exactly  analogous  way: 


£*(/)  = 


Mf) 
A0(f)  ' 


where  Ar(/)  is  the  amplitude  spectrum  of  the  reflected  Rayleigh  wave. 


(7b) 


Transmitted  Rayleigh  waves  are  recorded  by  keeping  source  and  receiver 
at  opposite  points  on  the  diameter  of  a  circle  centered  on  the  step  (see  Figure 
12).  Reflected  waves  are  observed  by  placing  source  and  receiver  at  equal 
angles  from  the  normal  to  the  strike  of  the  step.  Incidence  angle  is  measured 
from  the  strike  of  the  step. 


111.3  MODELING  RESULTS 
Vertical  Step  Model 

Figure  13  shows  seismograms  of  Rayleigh  waves  transmitted  across  a  3  mm 
vertical  step  at  incidence  angles  ranging  from  15°  through  90°  (normal 
incidence).  The  small  signal  arriving  about  1  psec  before  the  Rayleigh  wave  is 
an  S  wave  excited  at  the  step.  At  a  given  incidence  angle  the  waveforms  for  the 
case  of  upstep  and  downstep  transmission  are  essentially  identical.  Note  that 
not  only  do  peak  amplitudes  decrease  by  a  factor  of  two  to  four  or  more,  but 
also  that  the  waveform  changes  strongly  as  a  function  of  angle  of  incidence.  At 
normal  incidence  (0=90°)  transmission  coefficients  (Figure  14)  agree  well  with 
those  obtained  for  a  two-dimensional  step  using  theoretical  (numerical) 
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calculations  (Drake.  1972;  Martel  et  al.,  1977;  Fuyuki  and  Nakano,  19B4). 

At  normal  incidence  the  fraction  of  transmitted  energy  (Ey)  seldom 
exceeds  20%  within  the  usable  frequency  band.  This  represents  a  loss  of  over 
80%  of  the  incident  Rayleigh  wave  energy  to  reflected  surface  waves  and 
converted  body  waves. 

Reflected  Rayleigh  wave  seismograms  are  complicated  by  the  presence  of 
double  reflections,  separated  by  about  2  /nsec  (Figure  15).  In  Figure  15a.  the 
first  reflection  corresponds  to  the  Rayleigh  wave  reflected  from  the  bottom 
edge  of  the  step,  the  second  to  the  Rayleigh  wave  reflected  from  the  upper  edge 
after  propagating  up  the  face  of  the  step.  These  are  further  complicated  by 
multiple  scattering  and  Rayleigh  to  shear  scattering.  Reflections  from  a 
downstep  (Figure  15b)  are  much  larger  than  those  from  an  upstep  (Figure  15a). 
Also,  shapes  of  down-step  reflected  waveforms  are  different  from  those  of 
upstep.  They  do  not  exhibit  the  strong  incidence  angle  dependence.  Similar 
differences  between  upstep  and  downstep  scattering  of  Rayleigh  waves  have 
been  observed  in  two-dimensional  modeling  experiments  (e.g.  Nathman,  1980; 
Chapter  II). 

In  general  the  reflected  energy  is  very  small,  particularly  at  normal  or 
near-normal  incidence  (Figure  17).  At  most  angles  and  frequencies,  reflected 
energy  accounts  for  no  more  than  10-20%  of  the  incident  energy.  As  the 
incidence  angle  becomes  more  grazing  the  reflected  energy  increases  rapidly, 
especially  for  wavelengths  approximately  equal  to  three  times  the  step  height 
(/*0.3). 

It  is  clear  from  Figures  14  and  16  that  the  efficiency  of  body  wave 
scattering  is  strongly  dependent  on  the  angle  of  incidence.  In  the  3  mm  step 
model,  at  wavelengths  equal  to  one-half  the  step  height  (the  first  maximum  in 
Figure  14)  about  90%  of  the  incident  energy  is  scattered  to  body  waves.  As  the 
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incidence  angle  decreases  to  near  0°  (grazing  incidence),  this  figure  drops  to 
about  zero. 

Ramp  Model 

To  investigate  the  effect  of  a  more  gradual  change  of  elevation,  the  3  mm 
step  was  tapered  to  a  slope  of  60°  from  the  horizontal  (see  Figure  12). 

Distortion  of  both  transmitted  and  reflected  waveforms  is  strong,  particularly  at 
normal  and  near-normal  incidence  (Figure  17). 

Energy  transmission  coefficients  for  this  model  (Figure  18a)  are  generally 
higher  than  those  for  the  step  model  (Figure  14b)  at  frequencies  below  0.5. 
The  energy  minimum  near  /  «0.9  at  normal  incidence  is  considerably  smaller 
than  the  corresponding  minimum  for  the  step  model.  Reflected  energy  (Figure 
18b)  is  considerably  smaller  than  for  the  case  of  the  vertical  step  (Figure  14a). 
The  double  reflection  is,  however,  apparent  at  all  angles.  Interference  between 
the  two  reflections  and  the  scattered  body  waves  contributes  to  the  complexity 
of  the  seismograms. 

Note  that  in  the  case  of  the  ramp  model  the  frequencies  have  been 
normalized  to  a  step  height  of  h=3.0  mm;  in  fact  this  normalization  is  arbitrary. 
While  the  vertical  change  in  elevation  is  indeed  3.0  mm,  the  distance  from  the 
lower  edge  to  the  upper  edge  as  measured  along  the  ramp  surface  is  actually 

larger  than  this  (—r1 — ,  where  a  is  the  ramp  angle).  This  distance  may  also 
sina 

enter  into  the  scaling  relations.  Direct  comparisons  between  step  and  ramp 
transmission  and  reflection  coefficients  may  not,  therefore,  be  appropriate 
when  scaled  to  the  same  step  height. 

Azimuthal  Dependence  of  Scaling 

In  the  models  discussed  above,  the  character  of  the  transmitted  and 
reflected  waveforms  changes  most  rapidly  as  the  incidence  angle  decreases 


below  about  30°.  It  is  in  the  range  0°  to  30°  that  the  behavior  of  Rayleigh  wave 
propagation  departs  most  strikingly  from  that  predicted  by  two-dimensional 
theoretical  and  approximate  methods.  At  incidence  angles  greater  than  30° 
simple  empirical  relations  can  be  found  between  the  results  obtained  at  normal 
and  oblique  incidence. 

In  the  case  of  transmission  through  the  vertical  step  (Figure  14).  note  that 
the  energy  transmission  minimum  near  / =0.5  at  normal  incidence  moves  to 
higher  frequency  as  the  incidence  angle  decreases.  The  equivalent  effect  is 
observed  for  normally  incident  waves  when  the  step  height  is  decreased.  This 
suggests  that  the  effective  step  height  is  dependent  on  incidence  angle.  That  is, 
as  the  incidence  angle  decreases  the  incident  Rayleigh  wave  "sees”  a  shorter 
step.  The  functional  dependence  of  effective  step  height  on  incidence  angle  is 
estimated  by  finding  a  simple  functional  relation  between  incidence  angle  and  a 
particular  feature  in  the  transmission  curves.  The  normalized  frequency  at 
which  the  first  minimum  occurs  in  the  transmission  curves  is  approximated  well 
using  the  empirical  relation 


S  ,n\  _  /rrtn(90°)  / q „ 

/mm(9)  ~  sinQ  •  (  « 

where  /^(Q)  is  the  value  of  /  at  which  the  first  minimum  in  the  transmission 
curve  appears.  When  this  relation  is  applied  to  the  entire  transmission  curve 
the  curves  at  each  incidence  angle  can  be  made  to  coincide  ^.Figure  19).  In 
terms  of  the  effective  step  height  at  angle  0,  Equation  (8a)  implies 


h,ff(9)  =  hsinQ  (8b) 

To  a  first  approximation  the  transmission  coefficient  for  an  arbitrary  incidence 

angle  can  thus  be  estimated  from  the  normal  incidence  transmission  coefficient 

using  this  empirical  relation.  This  approximation  is  poor,  however,  at  incidence 


angles  shallower  than  about  30°. 
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Because  of  its  strong  incidence  angle  dependence,  the  reflection 
coefficient  for  the  step  is  difficult  to  model  with  a  simple  empirical  relation. 


111.4  IMPLICATIONS  FOR  THE  EARTH 

One  of  the  most  interesting  features  of  the  ultrasonic  results  is  that  both 
transmitted  and  reflected  Rayleigh  wave  pulse  shapes  are  severely  altered  as  a 
result  of  interaction  with  the  various  steps.  The  waveform  distortion  is  due  to 
the  combined  effects  of  frequency-dependent  energy  transmission  and 
reflection  and  interference  of  waves  multiply  reflected  within  the  step.  More 
complicated  input  signals,  such  as  dispers  wave  trains,  can  be  expected  to 
undergo  similar  kinds  of  angle-dependent  waveform  distortion  after 
propagating  across  a  step. 

The  effect  of  a  given  model  on  a  surface  wave  signal  is  determined  by 
convolving  the  desired  input  signal  with  the  impulse  response  of  the  model.  As 
noted  earlier,  the  Rayleigh  wave  propagating  across  the  aluminum  half-space 
represents  the  impulse  response  of  the  transducers  and  recording  system.  The 
impulse  response  of  a  given  model  for  waves  propagating  in  a  given  direction  is 
thus  the  result  of  deconvolving  the  transducer/recording  system  impulse 
response  fr  a  the  recorded  signal.  In  the  frequency  domain  we  may  write 

*.<•»>  ■  ¥rf  ■  <»> 

where  fm(0;«)  is  the  complex  transfer  function  for  the  desired  model  at 
incidence  angle  0,  S,n(9;u)  is  the  complex  Fourier  spectrum  of  a  Rayleigh  wave 
propagating  across  the  model  at  angle  0,  ( u )  is  the  half-space  Rayleigh  wave 

spectrum,  and  u  is  the  angular  frequency.  The  inverse  Fourier  transform  of 
this  quantity  yields  the  impulse  response  of  the  given  model. 
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As  a  test  of  the  effect  of  a  step  on  a  dispersed  wave  train,  a  chirped 
sinusoid  was  convolved  with  the  impulse  response  of  the  3  mm  vertical  step  and 
ramp  models  at  several  incidence  angles  (Figure  20).  At  some  incidence  angles 
and  frequencies  the  envelope  of  the  wave  train  drops  to  near  zero,  suggestive  of 
the  "beating"  phenomenon  frequently  observed  in  seismograms  of  Rayleigh 
waves  propagating  across  ocean  basins.  This  beating  is  usually  attributed  to 
multipathing  due  to  lateral  heterogeneities  along  the  propagation  path 
(Evernden  1953,1954).  At  least  part  of  the  amplitude  modulation  of  dispersed 
Rayleigh  wave  trains  observed  in  the  Earth  may  be  due  to  propagation  across 
sharp  structural  features  such  as  those  modeled  here. 

In  order  to  investigate  the  effect  of  simple  steps  on  Rayleigh  wave 
propagation  in  the  Earth,  the  physical  dimensions  of  the  model  are  scaled  to 
correspond  with  those  of  some  realistic  Earth  structure.  The  scaling  factor 
relating  the  physical  dimensions  of  the  Earth  to  those  of  the  model  is  given  by 
the  linear  relation  (White,  1965): 


7m  =*7.. 

where  the  scaling  factor  k  is  given  by 


k  =  (—)( —1 

and  7*.  Lx.  and  Vx  are  characteristic  times,  lengths,  and  velocities,  respectively 
in  either  the  model  (m)  or  the  Earth  (e ). 


.  In  general,  the  terrestrial  scale  length  appropriate  for  a  given  physical 
problem  is  considerably  greater  than  that  which  would  be  indicated  by  surface 
topography  alone.  In  discussions  of  terrestrial  scaling,  the  terrestrial  scale 
length  is  not  to  be  interpreted  literally  as  the  dimension  of  any  real  structure; 
it  is  simply  the  characteristic  dimension  of  a  highly  simplified  model.  The 
topographic  expression  of  the  physical  model  includes  the  combined  effects  of 
terrestrial  topography  and  deep  structure.  The  ultrasonic  step  models 
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described  in  the  first  section  of  this  chapter  represent  a  very  simplified 
approximation  to  complex  boundaries  of  vertical  heterogeneities  of  the  earth 
structure.  It  is  still  worthwhile,  however,  to  investigate  whether  models  produce 
effects  similar  to  those  observed  in  the  earth. 

We  take  as  the  input  signal  for  a  given  model  a  real  or  synthetic  terrestrial 
seismogram.  The  seismogram  is  resampled  at  a  rate  determined  by  the  model 
scaling  factor,  and  convolved  with  the  impulse  of  the  model.  The  result  of  the 
convolution  is  resampled  to  restore  the  seismogram  to  its  original  (terrestrial) 
time  scale. 

Because  the  usable  frequency  of  Rayleigh  waves  in  the  modeling 
experiments  is  limited  to  the  band  0.2— 2.0  MHz,  the  range  of  normalized 
frequencies  available  for  scale  modeling  is  also  limited.  The  range  of  step 
heights  and  wavelengths  which  may  be  modeled  reliably  given  the  band-limited 
nature  of  the  ultrasonic  pulse  is  shown  in  Figure  21. 

The  Ocean-  Continent  Margin 

The  incidence  angle  dependence  of  Rayleigh  wave  transmission  across  the 
ocean-continent  margin  can  be  seen  in  recordings  of  teleseismic  oceanic 
Rayleigh  waves  made  at  coastal  stations.  For  example,  Figure  22  shows  Rayleigh 
waves  recorded  at  WWSSN  stations  on  the  west  coast  of  the  United  States  from 
an  earthquake  in  the  Loyalty  Islands.  The  incidence  angle  is  taken  to  be  the 
angle  between  the  great-circle  path  and  the  continental  shelf.  For  paths  to 
BKS,  COR  and  LON,  the  angles  are  approximately  90°,  60°,  and  45°.  respectively. 
Among  the  continental  stations  in  this  example,  epicentral  distance  ranges 
from  67.2°  (BKS)  to  91.9°  (LON).  Great-circle  azimuths  vary  by  no  more  than 
about  1°  and  source  radiation  patterns  would  have  the  same  effect  at  these 
stations.  Differences  between  the  seismograms  recorded  at  the  continental 
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stations  are  therefore  primarily  due  to  variations  in  the  incidence  angle.  At 
normal  incidence  (BKS)  most  of  the  Rayleigh  wave  energy  is  confined  to  the  first 
dozen  or  so  cycles,  after  which  time  the  amplitude  decays  rapidly  to  a  small 
value.  As  the  incidence  angle  decreases  to  near-grazing,  the  initial  'packet'  of 
energy  broadens  and  the  rate  of  amplitude  decay  decreases.  This  trend  is 
consistent  with  that  observed  in  the  chirped  sine  synthetics  (Figure  20). 

Synthetic  Rayleigh  wave  seismograms  were  generated  for  a  pure  oceanic 
path  using  the  PEM-0  model  (Dziewonski  et  al.,  1975)  and  convolved  with  the 
inpulse  response  of  the  3  mm  ramp  models,  scaled  to  a  height  of  60  km  (Figure 
22b).  The  model-synthetic  seismograms  exhibit  azimuthal  behavior  similar  to 
that  of  the  real  seismograms  (Figure  22a). 

In  this  example  we  have  assumed  the  Rayleigh  waves  propagate  on  the 
great  circle  path  between  source  and  receiver.  In  fact,  the  velocity  contrast 
across  the  ocean- continent  margin  in  the  Earth  is  such  that  Rayleigh  waves  of  a 
given  period  traversing  the  margin  will  deviate  from  the  great-circle  path  by  an 
angle  governed  by  Snell's  law  (e.g.  Evernden,  1953;  Capon,  1970).  Using  the 
coast  and  ocean  models  given  by  Drake  and  Bolt  (I960)  we  find  that,  in  the 
period  range  of  interest  in  this  example  (15-40  sec),  Rayleigh  waves  incident  on 
the  coast  at  incidence  angles  greater  than  30°  are  refracted  laterally  by  no 
more  than  about  8°.  This  is  within  the  uncertainty  (±10°)  in  the  estimate  of 
incidence  angle  and  so  is  not  an  important  consideration  in  this  example. 

The  Tibetan  Plateau 

The  Tibetan  Plateau  is  often  regarded  as  an  isostaticaily-compensated 
block  of  uplifted  crustal  rock  with  an  average  elevation  of  5  km  (e.g.  Bird  and 
ToksSz,  1977).  According  to  the  ultrasonic  modeling  results,  we  expect  that 
Rayleigh  waves  crossing  the  margins  of  such  a  block  will  be  scattered  at  the 
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edges  of  the  block  and  the  waveforms  distorted.  To  investigate  the  effect  of  the 
Tibetan  Plateau  on  Rayleigh  waves  propagating  across  it  we  compare  real  and 
model-synthetic  Rayleigh  wave  seismograms  for  the  propagation  paths  shown  in 
Figure  23.  The  synthetic  seismograms  are  based  on  model  results.  The 
amplitude  and  phase  response  of  Rayleigh  waves  incident  at  the  upstep  and 
downstep  as  the  path  crosses  Tibet  are  taken  from  ultrasonic  model  results. 

The  Tibetan  Plateau  is  modeled  as  a  constant-elevation  feature  over  a 
uniform  half  space.  We  assume  Rayleigh  waves  follow  a  great  circle  path  from 
source  to  station.  We  determine  the  angle  the  great  circle  path  makes  with  the 
edges  of  the  plateau.  Then,  using  the  method  described  in  the  previous  section, 
model-synthetic  Rayleigh  wave  seismograms  are  generated  by  convolving  an 
input  waveform  with  the  impulse  response  of  the  step  model  for  each  of  the  two 
edges  (up-step  and  down-step)  with  the  appropriate  angles  of  incidence. 
Because  the  Tibetan  Plateau  is  wide  compared  to  a  Rayleigh  wavelength, 
interference  arising  from  multiple  reflections  from  opposite  edges  of  the 
Plateau  is  neglected.  Since  there  are  no  WWSSN  stations  on  the  northern  front 
of  the  Plateau,  we  use  the  Rayleigh  wave  recorded  at  KBL  as  an  input  signal  to 
the  lens  model.  The  regional  attenuation  model  of  Bird  (1976)  is  applied  to  the 
synthetic  data. 

The  wavelength  scaling  implies  that  the  Tibetan  Plateau  modeled  here  has 
an  "effective"  elevation  of  40  km,  which  is  higher  by  a  factor  of  eight  than  the 
true  elevation  of  the  plateau.  The  model  assumes  crust  and  mantle  are  uniform 
across  India  and  Tibet  and  the  "effective"  height  of  the  plateau  represents  the 
lumped  contribution  of  crustal  thickness  differences  as  well  as  topography.  In 
spite  of  this  simplification  the  agreement  between  real  and  modeled  waveforms 
is  generally  good,  particularly  at  longer  periods  (Figure  24).  The  observed 
(WWSSN)  seismograms  are  vertical  components  of  Rayleigh  waves  from  a 
presumed  explosion  at  Lop  Nor  (Bird  and  Toksbz.  1977).  Note  that  the  relative 


amplitudes  between  SHL  and  ND1  are  similar  for  both  the  observed  and  model- 
synthetic  seismograms.  No  attempt  has  been  made  here  to  model  the  small- 
scale  scattering  in  Tibet  and  the  Himalayas;  the  short  period  (<15-20  sec) 
components  of  the  waveforms  therefore  do  not  match  as  well. 

111.5  CONCLUSIONS 

Using  three-dimensional  ultrasonic  modeling  techniques,  transmission  and 
reflection  coefficients  as  a  function  of  wavelength  and  incidence  angle  are 
obtained  for  Rayleigh  waves  propagating  in  simple  step  models.  The  efficiency  of 
Rayleigh  wave  transmission  and  reflection  varies  strongly  with  both  incidence 
angle  and  wavelength.  At  some  incidence  angles  and  frequencies  as  much  as 
90%  of  the  incident  Rayleigh  wave  energy  is  scattered  into  body  waves. 
Transmission  coefficients  obtained  at  normal  (90°)  incidence  agree  with  those 
obtained  with  two  dimensional  modeling  methods  and  two-dimensional  finite- 
difference  methods.  An  empirical  relation  between  transmission  coefficients  at 
different  incidence  angles  is  obtained  for  the  case  of  a  vertical  step.  This 
relation  may  be  applied  to  two-dimensional  theoretical  (numerical)  calculations 
of  propagation  in  more  complicated  structures  to  obtain  estimates  of  Rayleigh 
wave  transmission  coefficients  at  oblique  incidence  angles. 

The  effect  of  simple  structures  on  Rayleigh  wave  propagation  in  the  Earth 
is  demonstrated  by  convolving  the  impulse  response  of  the  models  with  real  or 
synthetic  dispersed  wave  trains.  Some  observed  features  of  Rayleigh  waves 
crossing  an  ocean-continent  margin  and  the  Tibetan  Plateau  are  similar  to 
those  predicted  from  model  results.  The  method  may  be  applied  to  studies  in 
which  surface  waves  cross  regions  of  strong  structural  contrast.  In  Basin-and- 
Range  provinces,  for  example,  where  crustal  Rayleigh  waves  traverse  numerous 
steps,  strong  effects  can  be  expected  due  to  compounding  of  the  effects 
suffered  by  passage  through  a  single  step. 
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The  ultrasonic  model-synthetic  method  may  prove  useful  in  more  general 
problems  of  three-dimensional  surface  wave  scattering,  particularly  in 
situations  for  which  analytic  or  numerical  solutions  are  difficult  to  obtain. 
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FIGURE  CAPTIONS 

Figure  12.  Geometry  of  the  models  used  in  this  study.  Source  transducer  is 
located  at  point  S.  Reflected  or  transmitted  waves  are  recorded  at  points  R  or 
T.  Cross  sections  of  two  steps  used  are  shown  below.  The  source  and  receiver 
transducers  are  identical.  Each  transducer  is  made  of  a  piezoelectric  unit 
embedded  in  a  lucite  wedge  as  shown  at  the  bottom  of  the  figure.  This 
arrangement  maximizes  the  Rayleigh  wave  response.  The  wedges  are  coupled  to 
the  model  with  vacuum  grease. 

Figure  13.  Waveforms  of  Rayleigh  waves  transmitted  across  a  3mm  vertical  step. 
The  top  trace  is  the  input  signal,  obtained  by  putting  source  and  receiver  10  cm 
apart,  on  the  half  space.  The  small  pulse  arriving  1.0  /xsec  before  the  Rayleigh 
wave  is  an  S  wave  excited  at  the  step.  The  vertical  scales  on  the  seismograms 
have  been  expanded  by  the  scaling  factors  shown  at  the  start  of  each  trace.  A 
Transmission  in  the  up-step  direction.  B .  Transmission  in  the  down-step 
direction. 

Figure  14.  Energy  transmission  coefficients  for  the  Rayleigh  waves  in  Figure  13. 
The  horizontal  axis  is  normalized  frequency  (/  =  —),  where  h  is  the  step  height 

A 

and  X  is  the  wavelength.  Curves  are  shown  for  three  angles  of  incidence.  A 
Transmission  in  the  up-step  direction,  i?.  Transmission  in  the  down-step 
direction. 

Figure  15.  Waveforms  of  Rayleigh  waves  reflected  from  the  3mm  vertical  step. 
Note  the  secondary  reflection  corresponding  to  a  reflection  from  the  upper 
edge  of  the  step.  A  Reflection  from  the  upstep.  B:  Reflection  from  the 
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downstep. 

Figure  16.  Energy  reflection  coefficients  for  the  reflected  Rayleigh  waves  in 
Figure  15.  Note  the  change  of  vertical  scale  from  Figure  14.  A  Reflection  from 
the  upstep.  B:  Reflection  from  the  downstep. 

Figure  17.  Seismograms  of  Rayleigh  waves  transmitted  across  (A)  and  reflected 
from  (B)  the  3mm  60c  ramp.  Propagation  is  in  the  up-step  direction  The 
seismograms  in  (B)  were  high-pass  filtered  at  0.3  MHz  at  the  time  of  recording 
to  enhance  the  signal-to-noise  ratio. 

Figure  18.  Energy  transmission  and  reflection  coefficients  for  the  Rayleigh 
waves  in  Figure  17.  A:  Upstep  transmission.  B:  Upstep  reflection.  Note  the 
difference  in  the  amplitude  scales. 

Figure  19.  Energy  transmission  data  of  Figure  14a  after  application  of  the 
empirical  relation  of  Equation  (3).  The  transmission  curves  for  the  different 
incidence  angles  now  coincide. 

Figure  20.  The  effect  of  the  3  mm  vertical  step  (A)  and  ramp  (B)  models  on  a 
dispersed  wave  train.  The  input  signal  (top  trace)  is  a  linear  chirped  sinusoid, 
whose  frequency  increases  from  0.2  to  2.0  MHz  in  60  p.sec.  The  corresponding 
range  of  the  ratio  of  step  height/wavelength  (/ )  is  0.2  to  2. 

Figure  21.  The  region  in  h-X  space  in  which  the  ultrasonic  model  results  from 
the  3  mm  models  may  be  reliably  used.  Units  of  h  and  X  are  arbitrary.  In  the 
shaded  areas  repeatability  and  consistency  of  a  given  experiment  is  not 
assured.  The  usable  and  unusable  regions  are  demarcated  by  the  maximum  and 
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minimum  usable  normalized  frequencies.  The  box  shows  the  region  Allowed  for 
modeling  Rayleigh  waves  in  the  period  range  15-40  sec  (50-150  km  wavelength). 
Thus  step  heights  in  the  range  30-100  km  may  be  modeled  successfully  over 
this  period  range. 

Figure  22.  A.  WWSSN  vertical  component  recordings  of  the  Rayleigh  waves  from 
an  earthquake  in  the  Loyalty  Islands  (9/B/B1,  llh02m40.Bs;  21.4BB0  S.  169.600° 
E.  depth=31  km.  11^=5. 9.  Ms=6.2).  The  angle  at  which  the  great-circle  crosses 
the  continental  margin  is  shown  next  to  the  continental  stations.  B  Model 
synthetic  seismograms  using  the  3mm  60°  ramp  scaled  to  60  km  height. 

Figure  23.  Rayleigh  wave  paths  across  the  Tibetan  Plateau  from  an  explosion  at 
the  Lop  Nor  test  site  (14  October,  1970).  Reproduced  from  Bird  (1976).  Inset: 
Three-dimensional  model  of  the  Tibetan  Plateau  used  to  generate  model- 
synthetic  seismograms  for  the  paths  Lop  Nor-NDl  and  Lop  Nor-SHL. 

Figure  24.  Comparison  of  real  seismograms  (solid  lines)  with  model-synthetic 
seismograms  (dashed  lines)  using  the  lens  model  of  Figure  12.  All  seismograms 
have  been  low-pass  filtered  at  10  second  period. 
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